AD-A140  630 
UNCLASSIFIED 


MARS:  SUBSURFACE  PROPERTIES  FROM  OBSERVED  LONGITUDINAL 
VARIATION  OF  THE  3..IU)  AEROSPACE  CORP  EL  SEGUNDO  CA 
ELECTRONICS  RESEARCH  LAB  E  E  EPSTEIN  ET  AL .  15  MAR  84 
TR- 00841 4925-07 > -5  SD-TR-84-05  F/G  3/2 


Prepared  for 

SPACE  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 
Lot  Aafdet  Air  Force  Station 
P.O.  Box  92960,  Worldway  Porta!  Center 
Lot  Angeiet,  Calif.  90009 


V 

f  ' 

V" 


APRS  0S9A 


This  report  was  submitted  by  The  Aerospace  Corporation,  El  Segundo ,  CA 
90245,  under  Contract  Mo.  F04701-83-C-0084  with  the  Space  Division,  P.0.  Box 
92960,  Vorldway  Postal  Center,  Los  Angeles,  CA  90009.  It  was  reviewed  and 
approved  for  The  Aerospace  Corporation  by  D.  H.  Phillips,  Director,  Electronics 
Basearch  Laboratory.  Major  Robert  L.  Jones,  SD/7ASM,  was  the  project  officer 
for  the  Mission-Oriented  Investigation  and  Experimentation  (M0 IE)  Program. 

This  report  has  been  reviewed  by  the  Public  Affairs  Office  (PAS)  and  is 
releasable  to  the  National  Technical  Information  Service  (HTIS).  At  NTIS,  it 
will  be  available  to  the  general  public.  Including  foreign  nationals. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 
Publication  of  this  report  does  not  constitute  Air  Force  approval  of  the 
report's  findings  or  conclusions.  It  is  published  only  for  the  exchange  and 
stimulation  of  ideas. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  TNIS  PACE 


REPORT  DOCUMENTATION  PAGE 


in  4  1 4* 


SD-TR-84-05 


■i'-.i  E-  ■**  yt'.i’fZ 


wnsBf. 


*.  TITLE  MOM*; 

MARS:  SUBSURFACE  PROPERTIES  FRCM 
OBSERVED  LONGITUDINAL  VARIATION  OF 
THE  3.5-*ti  BRIGHTNESS  TEMPERATURE 


7.  autmoiv*) 

Eugono  E.  Epstein,  Bryan  H.  Andrew, 
Frank  H.  Briggs,  Bxuea  M.  Jskosky,  and 


(TTToNiF*  v>  i  rrrrr 


ANIMATION  NAM 

Tha  Aaroapaea  Corporation 
El  Sogundo,  Calif.  90245 


It-  CONTROLLINS  OFFICE  MAMS  AND  AODRKSS 

Space  Division 

Los  Angelas  Air  Force  Station 
Los  Angeles,  Calif.  90009 


>•  Tvn  OF  REPORT  A  PERIOD  COVERED 


•-  PERFORMING  ORS.  REPORT  RUMSER 


<<  ClIl,.?Tf \ 


IS.  REPORT  OATS 

15  March  1984 


Approved  for  public  release;  distribution  unlimited. 


17.  MSTRIRUTION  STATEMENT  (at  Ala  i 


tlmBHakH.  lit 


I*.  REV 


Mars,  Longitudinal  Brightness  T« 
Mars,  Millianter-Wave  Observatle 
Mars,  Subsurface  Properties 


iraturs  Variation 


s  at  3.5  as  of  the  disk-average  brightness  tenpnratnrn  of  Mars 
during  tbs  1978  opposition  cm  bn  reprossntod  by 


(Mars,  3.5  m ,  Jan/Fsb  1978)  -  INl*  10*eoe  (CML  ♦  JO^f*. 

("-J 


UNCLASSIFIED 


ggcuWTv  ci»wnc*Ttow  op  thu  e*es<os—  i 


»»•  Bv»58S7555SSr 


NT  A— TRACT  (SmSSmSST 


-(The  errors  cited  are  froa  the  internal  scatter;  the  eatlaated  absolute 
calibration  uncertainty  is  3X.)  This  longitudinal  variation  nust  be  taken 
Into  account  If  Mars  Is  to  be  used  as  a  calibration  source  at  allline ter 
wavelengths.  The  total  range  of  the  3.5^ea  variation  is  three  to  four  tines 
larger  than  both  the  2.8-ca  and  20-ya  variations.  This  unexpected  result 
can  possibly  be  explained  by  subsurface  scattering  froa  rocks  ^  1 .5-cn 
radius.  _  V 


PWCLASSIFIID 


PREFACE 


The  patience,  skill,  and  accuracy  of  Me.  Arlyn  Alonzo  are  greatly  appre¬ 
ciated.  The  work  of  E**  was  supported  by  NASA  Contract  No.  NASH  3255  and  Air 
Force  Contract  No.  F04701-83-C-0084  and  the  work  of  BHJ  by  the  Mars  Data 
Analysis  Program  of  NASA.  A  portion  of  the  research  reported  here  was  carried 
out  by  the  Jet  Propulsion  Laboratory,  California  Institute  of  Technology, 
under  contract  with  NASA.  The  telescope  operator  staff  at  the  NRAO  Kltt  Peak 
facility  are  thanked  for  their  able  assistance,  especially  during  the  good 
weather  when  the  data  acquisition  rate  became  quite  high.  Richard  Dltteon, 
D.  0.  Muhleman,  and  an  anonymous  referee  made  valuable  comments  on  the  manu¬ 
script.  We  thank  Stephen  J.  Kelhm  for  drawing  our  attention  to  the  Importance 
of  subsurface  scattering  and  for  useful  comments. 


Accession  For 

ImS  G"-.i  2 
r-TIC  TAB 

1 '  «A*  l *  1  Ill '  * 

Ju::Ur.  *  ‘ 


i 


PREFACE . 

I.  INTRODUCTION . 

II.  OBSERVATIONS . 

III.  RESULTS . 

IV.  DISCUSSION . 

V.  SUMMARY  AND  CONCLUSIONS 

REFERENCES . . . 


FIGURES 


1*  Disk-Average  3. 5-an  Martian  Brightness  Teaperatures 

vs.  Central  Meridian  Longitude.... .  12 

2.  Disk-Average  3 .5-no  Martian  Brightness  Teaperatures 

vs.  Central  Meridian  Longitude,  Bin  Averages . . .  13 

3.  Disk-Average  3 . 5-aa  Martian  Brightness  Teaperatures 

vs.  Central  Meridian  Longitude,  Adjusted  Data .  16 

4.  Disk-Average  3. 5-aa  Martian  Brightness  Teaperatures 
vs.  Central  Meridian  Longitude,  Adjusted  Data, 

Bin  Averages . . . . . .  17 

5.  Best  Fits  to  1978  Opposition  Martian  Distc-Average 
Brightness  Teaperatures  at  20  pa,  3.5  an,  and 

2.8  ca  vs.  Central  Meridian  Longitude . . .  19 

6.  Noontlae  Nadir-Viewing  Microwave  Brightness  Tenperature 
as  a  Function  o£  the  Surface  Theraal  Inertia  and 

the  Radio  Skin  Depth . . .  23 


5 
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INTRODUCTION 


Andrew,  Harvey  end  Briggs  (1977,1978)  aeasured,  and  Doherty,  Andrew,  and 
Briggs  (1979)  confined,  longitudinal  variations  in  Mars'  disk-average  2.8-cm 
radio  brlghtnees.  These  variations  indicate  variations  in  the  thermal  inertia 
and  radiometric  albedo  in  the  top  deelMter  or  so  of  the  Martian  soil.  We 
have  made  confirming  taeesurements  at  3.5  on.  We  use  the  3.5-mm,  2.8-ca, 
and  20-um  data  to  reach  conclusions  about  subaurface  properties. 

II.  OBSERVATIONS 


Observations  were  Mde  with  the  NRAO*  11-m  telescope  et  Kitt  Peak  during 
three  runs  in  1978:  January  9-12  and  29-31,  and  February  22-24.  The  observa- 
tione  on  January  9  were  obtained  during  an  antenna  pointing  calibration  run 
and  were  kindly  Mde  available  to  us  by  Dr.  B.  L.  Ullch.  Inclement  weather 
prevented  useful  measureMnts  on  January  11,  12,  and  31.  Table  1  describes  the 
Earth-Mars-Sun  geometry  at  the  tlMS  of  the  observations.  This  Mars  opposi¬ 
tion  was  quite  favorable;  Mars'  relatively  high -declination  permitted  observa¬ 
tions  through  most  of  the  night  and  Jupiter  and  Saturn,  which  served  as  refer¬ 
ence  sources,  were  also  at  relatively  high  declinations  and  were  conveniently 
spaced  -1-1/2  hours  in  R.A.  ahead  of  and  behind  Mars,  respectively. 

A  linearly  polarised  receiving  system  was  used;  the  double-sideband 
system  temperature  St  18.1  GHs  (3.5  m)  was  -350  K  and  the  IF  bandwidth  was 
5 00  MBs.  Standard  •on-on"  beaa  switching  with  a  6. 6-Ha  chopping  rate  was 


*The  National  Radio  Astronomy  Observatory  is  operated  by  Associated 
Universities,  Inc.,  under  contract  with  the  Retypes!  Science  Foundation. 


Table  1.  Parameters  of  the  Earth-Mars-Sun  Geometry  at  Times  of  Observations8 


Parameter 

Jan.  9-10 

Jan. 

29-30 

Feb. 

22-24 

Earth-Mars  distance  (AU) 

0.665  *  0.662 

0.664 

4 

0.668 

0.773 

♦  0.790 

Sun-Mars  distance  (AU) 

1.627  ♦  1.628 

1.642 

♦ 

1.643 

1.656 

♦  1.656 

Dg  (deg) 

14.0  ♦  13.8 

11.1 

♦ 

10.9 

9.2 

♦  9.2 

Ds  (deg) 

12.4  ♦  12.7 

15.6 

4 

15.8 

18.9 

♦  19.2 

Ag  +  180  (deg) 

219.9  ♦  219.4 

212.6 

4 

212.1 

206.7 

«•  206.4 

Ag  -  Ag  (deg) 

-11.4  ♦  -10. 2b 

4.7 

4 

5.9b 

21.3 

♦  22. 7b 

k 

0.991 

0.996 

0.961 

♦  0.956 

1 

11.2  ♦  10.0 

6.4 

4 

8.0 

22.8 

♦  24.1 

Central  Meridian 

Longitude  (deg) 

195  ♦  312 

347 

4 

138 

74 

♦  246 

*1)2  and  Dg  are  the  planetocentrlc  declinations  of  the  Earth  and  the  Sun,  respect¬ 
ively,  Ag  end  Ag  Are  the  planetocentrlc  right  ascensions  of  the  Earth  and  the  Sun;  k 
la  the  ratio  of  the  illuminated  disk  to  the  entire  disk;  and  1  is  the  angle  between 
the  Earth  and  the  Sun  as  seen  from  Mara. 


^Negative  and  positive  values  of  (Ag  -  Ag)  correspond  to  observations  with  the  sub¬ 
earth  point  in  the  afternoon  and  morning  hemispheres,  respectively. 
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employed.  The  been  throw  of  228  ere  eec  (about  2.6  half-power  beamwldths)  was 
In  azimuth  to  eliminate  the  effects  of  differential  atmospheric  emission.  An 
observation  was  2  or  3  minutes  in  duration,  depending  upon  whether  it  was 
taken  with  the  conventional  "FIVE"  or  "SEQUENCE"  technique.  (A  SEQUENCE 

observation  consists  of  a  sequence  of  on-on  pairs  of  readings  wherein  the 

telescope  la  directed  so  that  the  planet  is  measured  first  in  one  position  of 

the  switched  beam,  then  in  the  other  position.  A  FIVE  observation  consists  of 

short  SEQUENCE  observations  in  the  nominal  direction  of  the  planet  and  at  the 
four  cardinal  points  one-half  beamwldth  away;  it  takes  more  telescope  time, 
but  it  permits  the  telescope  pointing  errors  to  be  derived  and  compensated  for 
in  the  data  reduction.)  A  FIVE  observation  on  each  planet  was  always  taken 
first  to  adjust  the  telescope  pointing  for  the  subsequent  SEQUENCE  observa¬ 
tion. 

A  FIVE  and  a  SEQUENCE  observation  of  Jupiter  and  then  of  Mars  were  made 
through  the  night,  until  Jupiter  was  at  too  low  an  elevation,  at  which  dee 

Saturn  was  substituted.  Focus  checks  were  made  every  2  or  3  hours.  Atmo¬ 

spheric  attenuation  was  also  determined  at  similar  Intervals  via  "tipping". 
The  senith  optical  depth  ranged  from  0.04  to  0.10  throughout  the  three  observ¬ 
ing  runs,  save  for  a  few  hours  on  January  10  when  it  was  as  high  as  0.21.  The 

proximity  in  the  sky  of  the  reference  sources  reduced  the  deleterious  conse¬ 
quences  of  uncertainties  in  the  atmospheric  attenuation  corrections. 

Absolute  calibration  was  based  on  the  3.5-wm  brightness  temperature  of 
Jupiter  —  179.4  A  4.7  1  —  determined  by  Ullch  at  al.  (1980).  The  estimated 
overall  absolute  calibration  uncertainty  is  3Z.  The  changes  in  the  solid 
angles  of  the  planets  over  the  several  weeks  of  measurements  were  accounted 
for  la  the  data  reduction.  In  calculating  the  source  else  corrections  the 
antenna  pattern  was  taken  to  be  Gaussian  (though  not  constant  —  see  below) 


and  che  planetary  disks  were  assuaed  Co  be  uniformly  bright;  Saturn's  rings 
were  neglected.  However,  an  observed  Jupiter/Saturn  signal  strength  ratio  was 
detarmined  from  all  the  measurements  and  used  in  reducing  the  Mars  data  ob¬ 
tained  when  only  Saturn  was  available  as  a  reference.  Because  the  ring  in¬ 
clination  angle  changed  only  1?5  during  the  course  of  the  measurements  and 
because  the  3-mm  flux  from  Che  Saturn  system  changes  very  slowly  with  inclina¬ 
tion  angle  (Epstein  et  al.,  1984  and  Ullch,  1961),  we  neglected  any  change  in 
flux  from  che  rings.  Although  nonthermal  emission  because  of  radiation  belts 
produces  longitudinal  variations  of  Jupiter's  centimeter  and  decimeter  radia¬ 
tion,  this  nonthermal  ^emission  is  less  than  5Z  of  Jupiter's  total  emission 
shortward  of  1.5-cm  wavelength  (Berge  and  Guilds ,  1976);  the  rotational  vari¬ 
ability  of  the  nonthermal  emission  produces  an  even  smaller  percentage  varia¬ 
tion  in  the  total  esilsslon  shortward  of  l.S  cm.  Thus  at  3  mm  the  nonthermal 
emission  component  must  be  quite  small.  The  zonal  structure  of  the  visible 
disk  of  Jupiter  is  sufficiently  longitudinally  symmetric  that  we  do  not  expect 
significant  longitudinal  variations  of  Jupiter's  3-mm  thermal  emission. 
Similar  statements  apply  to  Saturn.  We  know  of  no  evidence  to  suggest  con¬ 
trary  expectations  at  3  mm. 

The  11-m  reflector  was  not  mechanically  stable;  this  fact  meant  that 
there  was  no  single  relationship  between  antenns  gain  and  antenna  attitude 
that  was  valid  at  all  times*  Similarly,  the  antenna  pattern,  as  manifested  in 
the  azimuth  and  elevation  half-power  beanwldths  determined  from  che  FIVE 
observations ,  was  not  constant  with  tins*  The  beamwidth  changes  from  night  to 
night  were  not  statistically  significant,  but  the  changes  between  che  three 
runs  were.  These  changing  beanwldths  ware  taken  into  account  when  subsequent¬ 
ly  reducing  the  FIVE  data.  To  process  the  data  we  did  assume  that  G(s),  the 
relation  between  gain  and  zenith  angle,  was  stable  for  half  nights  at  a  time; 
we  determined  G(s)  from  the  relevant  Jupiter  and/or  Saturn  data. 
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III.  RESULTS 


The  131  values  of  Che  disk-average  3.5-mm  Mars  brightness  temperatures 
are  plotted  against  Central  Meridian  Longitude  in  Fig.  1.  The  circled  symbols 
represent  data  for  which  Saturn,  not  Jupiter,  served  as  the  reference.  No 
attempt  has  been  made  to  estimate  the  uncertainties  of  the  individual 
values.  From  working  with  the  data  we  believe  that  the  major  source  of  scat¬ 
ter  is  the  antenna  gain  instability;  in  comparison,  uncertainties  In  inter¬ 
mediate  calibration,  source  size  correction,  and  atmospheric  attenuation 
correction  are  negligible. 

The  same  data,  but  grouped  into  longitude  Interval  bins  of  equal  width, 
are  shown  in  Fig.  2.  Although  the  scatter  is  large,  the  dependence  upon  CML 
Is  clear.  There  Is  an  unfortunate  paucity  of  data  near  the  indicated  peak 
(CML  -  320°).  A  one-harmonic  least-squares  fit  to  the  unblnned  data,  with 
equal  weights  assumed  for  all  131  points  (it  is  not  possible  to  assign  reli¬ 
able  weights  to  the  points),  yields 

TB(Mars,  3.5  mm,  Jan/Feb  1978)  -  198  +  10  cos(CML  +  42e)  K;  (1) 

±1  ±2  ±12 

the  standard  deviation  of  Che  residuals  is  8.6  K.  The  errors  cited  represent 
the  Internal  scatter  only;  the  estimated  absolute  calibration  uncertainty  is 
3X.  (Ha  explicitly  denote  the  epoch  of  observations  because  Mars'  brightness 
temperature  la  a  function  of  apoch  and  observational  aspect  (Cuzzi  and 
Muhlenan,  1972).]  A  two-harmonic  fit  yields 
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CENTRAL  MERIDIAN  LONGITUDE  (deg) 

Disk-Average  3.5-bo  Harden  Brightness  Temperatures  vs. 
Central  Meridian  Longitude.  The  circled  synbols  represent 
data  for  which  Saturn  served  as  the  reference;  Jupiter 
served  as  the  reference  for  all  other  data. 


TgOiars ,  3.5  on,  Jsn/Feb  1978)  - 

198  +  12  cos(CML  ♦  48*)  -  4  cos(2  CML  -  76*)  K;  (2) 

*1  ±2  *12  *3  *20 

because  the  standard  deviation  of  the  residuals  reduces  only  to  8.3  K ,  and 
because  the  saplltude  of  the  second  harmonic  term  Is  not  statistically  signi¬ 
ficant,  this  two-harmonic  fit  is  not  meaningful.  The  one-harmonic  fit  is 
shown  in  Fig.  2. 

These  results  indicate  that  if  Mars  is  to  be  used  as  a  calibration  source 
at  3  am,  the  CML  at  the  time  of  observation  must  be  taken  into  account  if 
precision  to  better  than  *>  *52  is  sought. 

Data  from  the  early  January  and  February  runs  overlap  in  the  interval 
195*  <_  CML  <_  246*;  in  this  interval  the  average  of  the  February  data  is  5.5  * 
0.82  lower  than  the  average  of  the  early  January  data.  Data  from  the  late 
January  and  February  runs  overlap  in  the  interval  72*  CML  140*;  the 
February  average  is  4.3  *  2.02  higher  than  the  late  January  average.  We 
attribute  these  discrepancies  to  the  combination  of  three  effects:  a)  the 
sub-earth  local  time  of  day  changed  by  about .  2  hours  over  the  three  runs , 
l.e.,  (Ag  -  Ag)  changed  by  ■  33*;  b)  the  polar  cap  was  rapidly  shrinking;  and 
c)  the  planetocentrlc  declination  of  the  Sun  (Dg)  increased  »  7*  as  the  de¬ 
clination  of  the  Barth  (Dg)  decreased  by  *  5*.  (Assuming  that  the  3 -mm 
brightness  temperatures  scale  with  heliocentric  distance  as  r”1^4  (Doherty  et 
al. .  1979)  means  that  the  brightness  temperature  decreased  only  0.42  during 
the  course  of  our  observations.]  Effect  (a)  means  that  the  observed  sub-Earth 
points  In  early  January  and  February  corresponded,  respectively,  to  early 
afternoon  and  to  late  morning  on  Mars.  Hence  we  expect  the  February  (late 
morning)  average  to  be  lower  than  the  early  January  (afternoon)  average. 
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Similarly,  «•  expect  the  February  average  to  be  lower  than  the  late  January 
average,  but  by  not  ae  much.  However,  the  decline  in  the  dlek-average  bright¬ 
ness  temperature  as  the  phase  angle  changes  is  likely  to  be  at  least  somewhat 
offset  as  the  cold  polar  cap  shrinks  and  the  Sun's  declination  lnereaaes. 
Predicting  numerical  values  for  the  combination  of  these  effects  would  require 
extensive  calculations,  and  the  values  would  be  model  dependent.  Instead,  we 
made  an  empirical  adjustment  for  these  affects  by  simply  scaling,  aeparately, 
all  of  the  early  January  data  and  all  of  the  late  January  data  to  bring  into 
agreement  the  respective  average  brightness  temperatures  with  the  February 

averages  in  the  respective  longitude  Intervals  of  overlap. 

% 

The  results  after  these  adjustments  are  shown  in  Figs.  3  and  4.  The 
overall  scatter  is  visibly  reduced.  The  one-harmonic  fit  to  the  unbinned  data 
(shown  in  Fig.  4)  is 


TgOfmrs,  3.5-nm  adjusted  data,  Jan/Fcb  1978)  »  198  +  10  cos(CML  ♦  10*)  K;  (3) 

*1  ±1  ±6 


the  standard  deviation  of  the  residuals  is  7.4  L  A  two-harmonic  fit  yields 

TgOfars,  3.5-mm  adjusted  data,  Jan/Feb  1978)  • 

198  ♦  11  cos(CML  +  12*)  +  2  cos(2  CML  ♦  29*)  K;  (4) 

±1  *2  *7  *2  *35 

the  standard  devl  n  of  the  residuals  is  7.2  K.  Again,  the  two-harmonic  fit 

is  not  meaning  u  one-harmonic  fit  to  the  adjusted  date  is  essentially 

the  seme  as  Bq.  s  >r  a  30*  phase  shift. 

For  comparison,  Doherty  at  al.  (1979)  obtained  the  following  two-harmonic 
fit  to  their  1978  opposition  2.8-cm  date: 
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(5) 


Tg(Mars,  2.8  cm,  Jan/Feb  1978)  * 

191.5  +  1.1  cos (CML  -  37*)  +  1.6  coa(2  CML  -  134*)  K. 

±0.4  ±0.8  ±38  ±0.8  ±18 

[Doharty  al.  corrected  chair  data  to  a  mean  heliocentric  dletance  of 
1.524  All  using  an  r"1^  scaling  factor;  we  have  reaoved  this  correction  froa 
Eq.  (5).  The  3.5-na  and  2.8-ca  data  ware  taken  at  alaost  identical  heliocen¬ 
tric  distances.]  The  higher  of  the  two  2.8-ca  peaks  occurs  at  CML  -  60*. 
l.a. ,  -  70*  after  the  3.5-aa  peak;  however,  the  3.5-aa  trough  and  the  deeper 
2.8-ca  ainlaua  both  occur  at  CML  -  160*  (see  Fig.  5).  The  eetlaated  absolute 
calibration  uncertainties  in  the  3.5-aa  and  2.8-ca  data  are  3Z  (6  K)  and  22 
(4K)  on  the  constant  teras  in  Eqs.  (1)  and  (3),  respectively;  therefore  the 
difference  between  the  asan  levels  of  the  3.5-ea  and  2.8-ca  data  is  not  signi¬ 
ficant.  There  is  a  difference,  however,  in  the  aaplltudes  of  the  curves, 
which  is  discussed  below. 

IV.  DISCUSSION 

The  Mars  rotational  curve  at  3.5  as  can  be  coapared  with  those  obtained 
at  other  radio  wavelengths  and  in  the  thermal  inf  fared  (Fig.  5).  The  Viking 
Or  biter  1  20-ya  observations  aede  at  the  saas  opposition  (Christensen  at  al.. 
1978;  P.  R.  Christensen,  personal  coeauni cation,  1981)  show  a  total  variation 
of  the  disk-average  brightness  teaperature  with  longitude  of  -  7  K  and  are 
represented  by 

Tb  (Mars,  20  im.  Doc  1977/Feb  1978)  - 
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Best  Vita  to  1978  Opposition  Martian  MsHlmi|i 
Brightness  Tonporatoros  at  20  vn,  3.5  an.  and 
2.8  ca  ts.  Central  Meridian  Longltnde.  The  20-pn 
result  (lq.  (8)]  la  froa  Christenson  et  al.  (1978), 
the  3.3-an  curve  (Bq.  (3)1  is  the  sane  as  in  Fig.  8, 
and  the  2.8-ca  result  (Bq.  (5)1  la  froa  Boherty 
et  al.  (1979). 


(6) 


239.4  -  1.7  cos(CML  -  67*)  -  2.4  cos(2CML  -  90*)  R; 

*0.1  *0.1  *6  *0.2  *3 

chc  standard  deviation  of  the  raalduala  ia  0.6  K.  Thaaa  data  refer  to  a  sub¬ 
spacecraft  local  time  of  ■  1/4  hour  before  noon. 

The  variations  correspond  to  predominantly  viewing  regions  of  high  ther¬ 
mal  inertia  and  low  albedo  or  regions  of  low  inertia  and  high  albedo,  with  the 
low-inertia  regions  coinciding  with  the  greatest  temperatures.  Low  thermal 
inertia  regions  are  characterised  by  a  large  amplitude,  shallowly  penetrating 
temperature  wave,  whereas  high,  thermal  inertia  regions  experience  e  lower 
amplitude  but  more  deeply  penetrating  temperature  wave.  This  effect  results 
from  the  control  of  the  surface  thermal  inertia  (defined  as  /KpC,  where  K  is 
the  thermal  conductivity,  p  is  the  density,  end  C  is  the  specific  heat  of  the 
surface  material)  on  the  surface  temperature,  such  that  a  region  of  low  ther¬ 
mal  inertia  will  have  higher  temperatures  near  noon  than  a  region  of  high 
thermal  inertia  (e.g.,  Ki offer  et  al.,  1977).  Because  Earth-based  observa¬ 
tions  are  generally  of  the  day  side  of  Mars,  these  near-noon  temperatures  are 
emphasised  (Christensen  et  al..  1978).  There  is  considerable  thermal  diver¬ 
sity  on  Mars  (Palluconl  and  Kief  far,  1981),  such  that  regions  of  low  or  high 
inertia  are  alternately  visible  as  the  planet  rotates,  producing  a  significant 
rotational  variation  of  brightness  temperature. 

The  rotational  curve  is  quite  different  at  radio  wavelengths.  Measure- 
manta  at  2.8  cm  (Andrew  at  al..  1977,  1978;  Doherty  et  al..  1979),  which  show 
a  peak-to-peak  variation  in  disk-average  temperature  with  amplitude  of  about 
5  K,  compared  to  about  20  K  at  3.3  am,  indicate  that  the  high  temperatures 
correspond  to  preferential  viewing  of  the  regions  of  high  inertia  (Jakosky  and 
Huhlenan,  1980).  This  difference  can  be  understood  in  terms  of  the  subsurface 
temperatures  and  the  penetration  depth  of  the  2.8-cm  radiation.  The  e-foldlng 
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electrical  skin  depth  Is  typically  ~10X ,  or  -30  ca  here  (see  Campbell  and 
Ulrlchs,  1969;  Muhleman,  1972;  or  Cuszi  and  Muhleaan,  1972),  such  that  the 
bulk  of  the  observed  energy  comes  from  beneath  the  diurnal  thermal  skin  depth 

of  5-10  cm.  Therefore,  the  observations  are  most  sensitive  to  the  temperature 

/ 

at  depth;  because  hlgh-lnertla  regions  have  a  greater  average  dally  tempera¬ 
ture  than  lov-lnertla  regions,  due  to  the  nonlinearity  of  the  Planck  function 
and  the  generally  lower  albedo  In  hlgh-lnertla  regions,  the  trend  Is  to  see 
higher  disk-average  temperatures  when  viewing  the  hlgh-lnertla  regions 
(Jakosky  and  Muhleaan,  1980).  The  dominant  effect  is  the  Planck  function 
nonlinearity.  To  maintain  radiative  equilibrium  over  an  entire  day,  there 
must  be  a  balance  between  the  lower  emission  rates  during  the  day  (because  of 
lower  daytime  temperatures  of  hlgh-lnertla  regions)  and  the  higher  emission 
rates  at  night  (caused  by  higher  nighttime  temperatures) .  Because  of  the 
nonlinearity,  the  Increase  In  the  nighttime  temperature  must  be  greater  than 
the  decrease  In  daytime  temperatures  to  produce  the  same  magnitude  of  effect 
on  the  emission  rates  (Jakosky  and  Muhleaan,  1980). 

It  la  because  the  Infrared  and  radio  observations  are  so  different  that 
the  3.5-mm  measurements  reported  here  are  of  Interest  (In  addition  to  obtain¬ 
ing  data  pertaining  to  Mars  as  a  millimeter-wave  calibration  source).  At  some 
wavelength  between  20  ye  and  2.8  cm  the  behavior  must  change  from  that  seen  at 
2.8  cm,  with  hlgh-lnertla  regions  having  higher  disk-average  temperatures,  to 
that  seen  at  20  ya ,  with  low-inertia  regions  having  higher  disk-average 
temperatures.  That  Is,  there  Is  likely  a  wavelength  where  the  disk-average 
brightness  temperature  does  not  vary  significantly  with  longitude.  The  wave¬ 
length  at  idiieh  this  occurs,  as  well  as  the  suite  of  rotational  curves  at 
nearby  wavelengths,  will  yield  inforaation  relevant  to  the  physical  nature  of 
the  surface  of  Mars,  as  will  be  discussed  below. 
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Although  properly  calculating  the  rotational  curve  for  Mars  at  radio  and 
Infrared  wavelengths  Is  a  formidable  teak,  even  for  the  simplest  cases,  some 
simple  calculations  can  be  done  to  identify  the  wavelength  at  which  the  cross¬ 
over  in  behavior  nay  occur.  Figure  6  shows  contour  plots  of  the  modelled 
noontime  nadir-viewing  microwave  brightness  temperature  as  a  function  of  the 
thermal  Inertia  and  the  radiation  penetration  e-foldlng  skin  depth  for  several 
sets  of  assumptions  regarding  amlsslvitles  at  radio  and  Infrared  wave¬ 
lengths.  These  results  were  obtained  by  calculating  the  diurnal  surface 
temperature  variation  for  a  given  set  of  thermal  properties,  extrapolating  the 
temperature  variation  to  depth,  and  integrating  the  emission  of  energy  over 
depth  at  local  noon,  using  a  technique  similar  to  that  described  by  Muhleman 
(1972)  and  Jakosky  and  Muhleman  (1980).  Each  calculation  is  done  assuming  a 
homogeneous  surface,  with  properties  uniform  with  depth;  thermal  inertia  and 
surface  albedo  were  taken  to  be  related  as  described  by  Kleffer  et  al. 
(1977).  In  Fig.  6a,  the  infrared  ealssivlty  eIR  was  taken  to  be  1.0,  and  the 
microwave  emlsslvlty  e^  waa  a  constant;  In  Fig.  6b,  eIR  was  taken  to  be  cor¬ 
related  with  albedo,  and  hence  Inversely  with  inertia,  as  per  the  results  of 
Christensen  (1982);  in  Fig.  6c,  ejR  correlates  with  albedo  and  a^  was  taken  to 
be  related  to  Inertia  as  indicated  by  the  correlation  of  radar  erosa-section 
and  inertia  (Jakosky  and  Muhleman,  1981).  These  assumptions  represent  the 
range  of  possible  variations  of  surface  properties,  exclusive  of  roughness  and 
scattering  effects, which  are  discussed  further  below. 

The  penetration  depth  at  which  the  noontime  temperature  varies  the  least 
with  thermal  inertia  (and  hence  with  longitude  as  the  planet  rotates)  is  that 
for  which  the  rotational  curve  would  be  expected  to  show  the  least  variation 
of  temperature  with  longitude.  From  Fig.  6,  this  situation  occurs  for  a 
penetration  depth  of  between  0.6  and  1.6  cm,  depending  on  the  assumptions 
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THERMAL  INERTIA 


RADIO  SKIN  DEPTH  (cm) 


rif.  6.  Hoontlat  Madlr-F  lowing  Nietcmw  Brightness  Ta^ntun 
as  a  Function  of  tho  Surface  Thermal  Inertia  and  the 
tadlo  Skin  Depth  (In  cn).  (a)  As tuning  constant  infrared  end 
alerowave  ealsslvlty.  (b)  Attuning  Infrared  ealsslvlty 
varies  with  albedo  (and,  hence,  with  Inertia). 

(e)  Sane  at  (b),  but  also  allowing  alcrowave  ealsslvlty 
to  be  correlated  with  thernal  inertia.  See  test  for  details. 


nade.  In  Che  absence  of  roughness  and  scattering  effects,  the  wave¬ 
length  X  and  the  penetration  skin  depth  6  are  related  by  6  ■  nX,  where  n  Is 
between  about  5  and  20  for  particulate  surfaces  (Campbell  and  Ulrlchs, 
1969).  Thus  the  cross-over  point  should  occur  for  an  observing  wavelength 
between  about  0.2  and  3.2  mm.  For  material  properties  like  that  of  the 
Moon  (n  *  13;  Muhleman,  1972;  Gary  and  Kalha,  1978).  the  corresponding  wave¬ 
length  range  would  be  0.3  to  1.1  an. 

Surface  roughness  may  affect  these  results  in  two  ways.  First,  centi¬ 
meter-sized  and  larger  surface  rocks  may  perturb  the  situation.  If  rocks  are 

er 

present  In  large  enough  numbers  to  affect  the  bulk  surface  thermal  inertia,  as 
Indicated  by  the  radar/thermal  analysis  of  Jakosky  and  Muhleman  (1981),  then 
the  above  calculations  may  not  be  adequate;  they  are  sufficient  If  the  rock 
abundance  were  to  not  vary  significantly  from  place  to  place  and  If  the  ther¬ 
mal  Inertia  were  determined  by  the  properties  of  the  fine  component  of  the 
surface,  as  suggested  by  Christensen  (1982).  In  the  former  case,  the  observed 
flux  will  be  Che  sum  of  that  from  the  fine  material  and  that  from  the  rock, 
weighted  by  the  fraction  of  the  surface  covered  by  each.  The  penetration 
depth  for  which  the  rock  and  fine  components  of  the  surface  have  the  same 
brightness  temperatures  at  noon  corresponds  to  the  wavelength  for  which  rota¬ 
tional  variation  of  disk-average  temperature  will  be  small  or  negligible;  for 
flat  rocka  larger  than  the  thermal  skin  depth,  under  assumptions  equivalent  to 
Chose  in  Fig.  6a,  this  squall ty  occurs  for  a  penetration  skin  depth  of  about 
1.3  cm,  within  the  range  discussed  above. 

Second,  surface  roughness  on  the  scale  of  the  observing  wavelength  or 
subsurface  scattering  will  modify  the  discussion  of  skin  depths*  In  the 
surface  roughness  case  •  observed  energy  will  have  been  scattered  at  the  sur¬ 
face,  effectively  having  been  emitted  from  shallower  regions  than  expected 
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based  purely  on  che  material's  electrical  properties.  The  effective  value  for 
n  will  then  be  smaller  than  the  range  of  5  to  20  discussed  above.  There  is  no 
direct  information  on  the  structure  of  the  Mars  surface  on  the  scale  of  inter¬ 
est  here,  about  0.2  to  3  mm,  except  at  the  Viking  landing  sites  where  no 

significant  roughness  on  these  scales  is  seen  (e.g. ,  Mutch  et  al.,  1976 
a,b).  Scattering  by  subsurface  rocks  will  have  che  same  effect,  such  that 
emission  will  be  from  shallower  regions  and  che  effective  value  for  n  will  be 
smaller.  Kelhm  (1982)  calculates  for  che  lunar  case  that  reasonable  sub¬ 
surface  rock  distributions  can  decrease  n  by  more  chan  a  factor  of  two  over 
the  homogeneous  case.  Because  of  the  subsurface  temperature  structure  near 
noon,  che  effect  of  decreasing  n  Is  to  increase  the  brightness  temperature. 
This  point  is  discussed  further  below. 

In  light  of  the  above  discussion,  the  cross-over  from  the  behavior  of  the 

rotational  curve  seen  at  radio  wavelengths  to  that  seen  In  the  infrared  should 

occur  at  a  wavelength  within  che  range  of  about  0.2  to  3  mm,  assuming  a  homo¬ 
geneous  surface  and  subsurface.  The  wavelength  at  which  it  occurs  will  yield 
information  on  the  Mars  electrical  properties  and  roughness-related  emisslvlcy 
effects.  The  3.5-mm  observations  presented  here  show  behavior  more  like  that 
seen  at  2.8  cm  rather  than  that  seen  at  20  urn.  The  lower  disk-average  temp¬ 
eratures  coincide  with  the  low-inertia  regions  near  120°  longitude  (the 
Tharsls  region),  while  the  higher  temperatures  generally  correspond  to  the 
high-inertia  regions  over  most  of  the  remaining  longitudes.  The  uncertainty 
In  the  measurements  aey  account  for  not  being  able  to  Identify  in  the  3.5-mm 
data  all  of  the  longitudinal  structure  seen  in  the  thermal  Inertia  map  of 
Palluconl  end  Kieffer  (1981). 

The  difference  between  the  absolute  levels  of  the  3.5-mm  and  2.8-cm 
curves  Is  not  statistically  significant  (see  Sec.  III).  The  similarity  in 
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levels  results  from  the  fact  that  both  wavelengths  are  sensitive  to  eaisslon 
from  depths  comparable  to  or  greater  than  the  thermal  diurnal  shin  depth  (on 
average,  about  4  to  6  cm;  Kieffer,  1976),  and  again  indicates  that  3.5  am  is 
longward  of  the  cross-over  point. 

The  difference  between  the  amplitudes  of  the  3.5-mm  and  the  2.8-cm  curves 
is  most  important.  For  a  homogeneous  subsurface  one  would  expect  the  ampli¬ 
tude  of  the  rotational  curve  to  generally  decrease  with  decreasing  radio 
wavelength  until  the  cross-over  point  is  reached,  after  which  the  amplitude 
would  increase  (but  out  of  phase  with  the  radio  curve)  with  further  decreasing 
wavelength.  That  the  amplitude  increases  in  going  from  2.8  cm  to  3.5  mm 
indicates  that  subsurface  scattering  may  be  important.  That  the  high  3.5-om 
brightness  teoiperatures  generally  coincide  with  the  high  temperatures  of  the 
2.8-cm  curve  and  with  the  high-inertia  regions  indicates  that  the  subsurface 
scatterera  are  located  in  the  hlgh-lnterla  regions.  To  effectively  scatter 
the  3.5-mm  but  not  the  2.8-cm  waves,  the  scattering  efficiency  must  be  larger 
at  the  shorter  wavelength;  the  scattering  results  presented  by  Hansen  and 
Travis  (1974)  indicate  this  to  be  the  case  if  the  scatterers  are  smaller  than 
about  1.5— cm  radius.  The  distribution  of  surface  rocks  seen  at  the  Viking 
lander  sites,  however,  peaks  near  5  -  10  cm  radius  (e.g..  Hutch  et  al . ,  1976a, 
b;  1977),  and  has  a  paucity  of  centimeter-sized  rocks  (Mutch  et  al.,  1977). 
It  is  conceivable  that  centimeter-sized  roughness  elements  on  the  surfaces  of 
the  larger  rocks  are  responsible  for  the  scattering,  or  that  the  lander  sites 
are  not  typical  of  the  entire  planet;  that  this  latter  idee  is  plausible  is 
suggested  by  the  differences  between  the  Viking  Lander  1  region  and  other 
locations  in  thermal  and  radar  properties,  both  of  which  are  also  related  to 
rock  abundance  (Jakosky  and  Muhleman,  1981). 
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The  requirement  chat  che  high-inertia  regions  contain  additional  scatter- 
ers  is  consistent  with  the  hypothesis  that  the  thermal  inertia  of  the  surface 
is  largely  controlled  by  the  abundance  of  surface  rocks  (Jakosky  and  Muhleman, 
1981),  but  Inconsistent  with  the  relatively  constant  rock  abundances  Inferred 
by  Christensen  (1982).  Christensen's  Inference  is  model -dependent,  however, 
in  that  the  results  obtained  would  differ,  depending  on  the  size  of  rock 
assumed  to  be  present;  the  correlation  of  rock  abundance  with  inertia  that  he 
presents  may,  in  fact,  indicate  this  to  be  the  case. 

V.  SUMMARY  AND  CONCLUSIONS 

The  extensive  3.5-mm  measurements  reported  here  show  a  variation  in  Mars' 
brightness  temperature  with  Central  Meridian  Longitude  that  is  generally  in 
phase  with  the  variation  at  2.8  cm,  and  opposite  in  sign  from  the  variations 
at  20  um.  This  phase  result  is  not  unexpected  because  3.5  mm  is  longer  than 
the  wavelength  at  which  the  phase  behavior  is  expected  to  change. 

Unexpected,  however,  is  the  result  that  the  3.5-mm  rotation  curve  ampli¬ 
tude  is  larger  than  the  amplitudes  at  both  20  pm  and  2.8  cm.  This  result  can 
be  explained  by  subsurface  scattering  from  rocks  smaller  than  1.5-cm  radius. 
A  correlation  of  subsurface  scatterers  with  the  location  of  the  high-thermal 
inertia  regions  would  be  consistent  with  the  hypothesis  that,  predominantly, 
the  rock  abundance  determines  the  thermal  Inertia. 

Because  the  amplitude  of  the  3.5-mm  rotation  curve  is  larger  than 
expected,  and  because  the  3.5-em  measurements  were  hindered  by  Instrumental 
difficulties,  further  millimeter-wave  measurements  are  important  —  at  •  3  an 
to  confirm  the  present  results  and  at  <  1  an  to  better  understand  the  nature 
of  Mars'  subsurface.  Additionally,  more  detailed  radiative  transfer  models 
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are  justified  to  better  account  quantitatively  for  the  present  behavior  and 
thereby  constrain  the  scattering  models  and  to  permit  adjustment  for  the 
seasonal  effects  discussed  in  Section  III.  Such  observations  and  models  are 
now  being  planned. 
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